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ABSTRACT 


Emlssl vities were calculated from the Nimbus 5 Electrically 
Scanning Microwave Radiometer (ESMR) over 25 km grid cells for the 
period September 1973 through May 1975 for the southern Great Plains 
including the western two-thirds of Kansas and Oklahoma and northwest 
Texas. These emissivities, normalized for seasonal temperature 
changes, were in excellent agreement with theory and measurements made 
from aircraft and truck sensors at the 1.55 cm wavelength of ESMR. 
These emissivities were related to crop moisture conditions of the 
winter wheat in the major wheat producing counties of the three 
states. High correlations were noted between emissivity and an ante- 
cedent precipitation index (API) used to infer, soil moisture for 
periods when the soils were essentially bare. The emissivities from 
ESMR were related through API and actual crop condition reports to 
progress of fall planting, adequacy of crop moisture for stand 
establishment, and periods of excessive moisture that necessitated 
replanting. Periods of prolonged frozen soil in the winter were 
observable at several grid points. The average emissivities of the 
canopy/soil surface during the maximum canopy development times in the 
spring showed a good agreement with moisture stress inferred from 
rainfall and yield data. Discriminant analyses of the emissivities 
for both rainfall and API produced probabalistic relationships of 
total rainfall and maximum rainfall for given sequences of spacecraft 
for early warning of crop conditions is strongly supported by the 
research. 
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Introduction 

Large scale crop condition monitoring is severely hampered by the 
lack of information on soil moisture conditions. Conventional soil 
moisture measurement techniques are simply too unwieldy and time- 
consuming to provide the density of data required. Even in the Great 
Plains, modeling of soil moisture from meteorological, geophysical, 
and crop information is not practical, due to the very limited input 
data available on a real-time basis. Yet, moisture information is 
vital in any effort to monitor crops such as winter wheat over large 
areas for conditions that impact cultural practices, growing 
conditions and, ultimately yield. 

Passive microwave remote sensors have the capability to provide 
useful crop moisture information with sufficient time and space 
resolution. Several studies have demonstrated that brightness temper- 
atures from sensors such as the Electrically Scanning Microwave Radio- 
meter (ESMR) and the Scanning Multi frequency Microwave Radiometer 
(SMMR) are highly correlated with significant rainfall events over 
large agricultural areas. The resolution is of the order of 25 km 
spatially and every two or three days temporally, which is entirely 
adequate for large scale crop moisture monitoring. This proposed 
research is to develop an early warning screening program for moisture 
deficiencies or excesses at planting time and critical growth periods 
and an index of accumulated crop moisture and plant stress for the 
wheat areas of the Great Plains. 
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Passive Microwave Remote Sensing 

Investigations by Cl hi ar and Ulaby (1975), Schmugge, et aj_. 
(1974, 1976a, 1976b, and 1977), Schmugge (1976 and 1977), and Newton 
(1977) demonstrated that the emitted radiation at microwave frequen- 
cies is a function of the moisture content of the emitting soil 
layer. Basically, air and dry soil have a very low dielectric con- 
stant, while that of water is the highest of naturally occurring abun- 
dant substances. As water is added to the soil, the dielectric cons- 

tant of the soil, air, and water mixture increases, with a resulting 
decrease in the emissivity. The emissivlty is related to the radia- 
tion received at the sensor antenna through the simplified relation- 
ship: 

T U T EL (D 

where T g is the radiation received, also termed the brightness 
temperature since it is linearly and directly proportional to the 
actual temperature of the emitting layer T The emissivity is **. 
At wavelengths of 1.55 cm, the dry soil emissivity will normally be in 
the .92 to .95 range, while emissivities in areas receiving heavy 

rains will be as low as .74. At these relatively short wavelengths, 

the emitted radiation is absorbed and reflected by surface roughness 
and vegetative cover, which vary considerably from one area to the 
next. Attempts to quantitatively map soil moisture with 1.55 cm pas- 
sive microwave remote sensors have not been particularly encouraging 
(for example, see Meneely, 1977). However, time series of brightness 


temperature for the same sensor footprint areas (wiis^e surface rough- 
ness and vegetative cover variations are significantly reduced) show a 
high correlation with rainfall history and infrared soil moisture for 
essentially bare soils (see McFarland and Blanchard, 1977; Theis, 
1979; Blanchard, 1981a; and Blanchard, 1981b). 

In these studies in the Great Plains winter wheat areas, the 
brightness temperatures from the Electrically Scanning Microwave 
Radiometer (ESMR) on the Nimbus 5 spacecraft, and the rainfall and 
temperature records from the Climatological Data were objectively 
analyzed to a 25 km grid using a modified Barnes exponential weighting 
function. The grid established for the intensive study area of 
Oklahoma is showr) in Figure 1. This grid is based on a polar 
stereographic map projection, true at 35°N, in order to accept 
latitude/longitude coordinates for input data. With this grid and the 
objectively analyzed values at each grid point, problems with missing 
data, variable spatial and temporal densities of input data, and data 
management are simplified. 

Investigations with ESMR 

Twenty-seven grid locations were selected for analysis of the 
ESMR data. The location of these grid cells in Kansas, Oklahoma, and 
Texas is shown in Figure 1. Each grid cell represents a 25 km square 
area in a county with a substantial acreage of winter wheat. The 

wheat acreage in each ot the approximately equal area counties for the 
winter wheat in 1973-74 and 1974-75 is in Table 1. 

The ESMR data set consisted of brightness temperature obser- 
vations from September 5, 1973 to May 30, 1975. Approximately 260 
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Figure ],. Location of Grid Cells for ESMR Analysis. 
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Table 1. Winter Wheat Acreage and Yield in the ESMR Grid Cell 


Counties. 

Kansas 


District 

Grid 

. County 

Northwest . 

102 

Sherman 

177 

Thomas 

West Central 

106 

Greeley 

s'. 

204 

Gove 

. 

256 

Ness 

Southwest 

133 

259 

Kearney 

Ford 

North Central 

352 

• Osborne 

# 

453 

Ottawa • 

Central 

356 

Barton 


456 

McPherson 

South Central 

361 

Barber 


410 

, Kingman 


459 

Sedgwick 

, 

• 

, ' Texas 

Northern High 


Lipscomb 

Plains 

240 

Northern Low 

* 

* 

Plains 

299 

Hardeman 



Oklahoma 

Northwest 

213 ’ 

Beaver’ 


288 

Harper 

. 

291- 

.. Ellis 

_West Central 

343 

Custer 

Southwest - 

296 

Greer 


372 

Kiowa 


424 

Cotton 

North Centr.al 

413 

• Alfalfa 

• 

515 

Noble 

Central 

442 

Kingfisher 

• 471 

Grady 


(to convert to bu/ac, divide by 0.0673) 


Planted Acreage Yield 

(1000s ha) (T/ha) 

1973-7A 1974-75 1973-74 1974-75 


79 

79 

2.4 

2.0 

94 

93 

2.3 

2.1 

•78 

83 

2.3 

1.6 

57 

58 

2.1 

2.5 

75 

80 

. 1.8 

2.0 

58 

60 

2.2 

• 1.6 

99 

99 

1.7 

1.9 

58 

. 60 

1.4 

• 1.9 

63 

68 

2.0 

2.0 

95 ’ 

99 

1.3 

• 1.9 

94 

. 103 

1.9 

2.0 

69 

72 

1.7 

1.8 

96 

101 

1.8 

1.9 

106 

113 

1.8 

1.8 


49 

60 

°»&»7 

1.2 

52 

64 

0.9<6 

1.5 

122 

131 

0.6 

0.9 

68 

71 . 

1 .4 

1.2 

48 

51 ‘ 

0.9 

1.0 

100 

. 113 

1.7 . . 

1.9 

40 

42 

1.1 

1.5 

107 

118 

1.3 

1.7 

64 

66 

1.7 

1.4 

125 

• 133 

i.9 

1.9 

60 

58 

1.7 

1.9 

115 

126 

1.4 

1.7 

39 . 

41 

1.5 

1.9 
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observations of ESMR brightness temperatures were Included in the 
analysis for each grid cell. The only significant break In the period 
of record occurred between June 8, 1974 and August 20, 1974 wft*n only 
two days with ESMR coverage were noted In the central portion of the 
grid. For the period September 5, 1,973 to May 30, 1975 ESMR bright- 
ness temperatures were available, on the average, every 2.4 days. 

The analysis grid is the same 25 x 25 km grid used In earlier 
phases of this study. Daily maximum and mi nlmum air temperatures, 
daily precipitation, and snow depth from the NOAA Climatological Data 
were objectively analyzed to each grid location. A modified Barnes 
exponential weighting function using the seven nearest observations 
was used for the objective analysis. The ESMR brightness temperatures 
from magnetic tape provided by the NASA Goddard Space Flight Center 
were also objectively analyzed to each grid with the same function. 
Only those observations within 35 degrees from nadir were used, which 
corresponded to a resolution ranging from 25 km at nadir to about 50 
km at 3b degrees scanning angle. Data were used when the ESMR 
coverage occupied only a portion of this grid; not all grids have ESMR 
values for any given day. 

The emissivity model used in the investigation is: 

' * V T A 

where e is the normalized emissivity 

Tg is the brightness temperature from ESMR 

T a is the daily maximum air temperature at the grid 

The emissivity obtained is an approximation of the true 
emissivity layer; a necessary approximation in order to remove the 
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effects of daily and seasonal temperature changes. The resulting 
emissivities are, however, very near the expected values and the 
seasonal temperature trend is effectively removed from the brightness 
temperatures, as shown in Tables 2 through 4. 

Emissivities at 1.55 cm for smooth surfaces are (Schmugge et al . 

1977) 


Water at 

20° 

0.40 

0.40 

Dry soil 



0.94 

Wet soil 

above 

field capacity 

0.60 

Pure ice 



0.92 


The emitting surfaces for ESMR observation however are rough, 
which increases the emissivity. Choudhury et a I. (1979) reported the 
effect of roughness as: 

■ (1- )(l-exp" h ) 

where a is the change in emissivity from the smooth surface 
emissivity, , for a roughness h. Their data showed .the best 

correlations for observed emissivities over rough fields at the 
Phoenix site (aircraft observations) were from a roughness of 0.6, If 
this figure is used, the emissivity for a moist, rough terrain is 
0.74, which coincides well with the lowest emissivities observed in 
the ESMR data sets. 

Correlations of Emissivity with Antecedent Precipitation 

To correlate the normalized emissivity with the rainfall history 
of each test grid cell, an antecedent precipitation index (API) was 
used (McFarland, 1975; McFarland and Blanchard, 1977; Blanchard, et 
aj_. 1981). The API model used is: 

API 1 = AP I i _i k(t) + (r i ) 0 - 891 (1) 
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Table 2. Fall Correlations of API and Emissivity. 


GRID 

Mean 

Emissivity 

Mean 

API 

Grid Cell 
R-Squared 

Averaged 
Grid Cell 
R-Squared 

102 • 

.92 

’ .40 

.29 

.30 

106 

.91 

.43 

.46 

' .47 

133 . 

.92 

.29 

.17 

.21 

177 

.92 

.39 

.30 

. .31 

204 

. .92 

'.45 

.43 

.42 

213 

.91 

.61 

.27 

.29 

240 

.92 • 

.78 

.28 

.28 

256 

.90 

.47 

.54 

.50 

259 

.90 

.64 • 

.49 . 

.51 

288 

.90 

• 1.34 

.45 

• ..46. 

291 

.92 

.70 

.32 

,34 

' 296 

.89 

• .95 

.49 

• .53' 

299 

.89 

.95 

. .49 

.54 ‘ 

343. 

.89 

.87 

.57 

.56 

352 

. .90 ‘ 

.85 

.37 

.40 

356 

.87 

• 1.03 

.47 

.48 

361 

.90 

1.12 

.54 

.56 

372 

.90 

.93 • 

.31 ‘ 

.35. 

410 

.88 ' 

1.17 

.44 

.44 

413 

.88 

1.12 

.57 

.60 

424 

.90 - 

.98 • 

.52 

.40 

442 

.89 

1.07 

.38 

.41 

453 

.89 

1.16 

.31 

'.34 

456 

' .88 

1.20 

.44 

.43 

459 

.88 

1.17 

;41 

.44 

471 

.90 

.93 

.19 

.18. . 

515 

.90 

1.20 

.17 

.21 
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Table 3. Summer Correlations of API and Emissivity. 


GRID 

Mean 

Emissivity . 

Mean 

'API 

R 2 

102 

.94 

.65 

.37 

106 

.92 

.47. 

.55 

133 

.92 

.79 

.59 

177 

.93 

.52 

.45 

204 

.93 

.77 

.54 

213 

.92 

.56 

.36 

240 

.92 

.58 

.29 

256 

.92 

.59 

.42 

259 

.92 

1.13 

.28 

288 

.91 

' .94 ; 


291 

, .91 

.85 

. - 16 

296 

.90 

1.13 * 

.39 

' 299 

• .90 

1.10 

.42 

343 

.90 

1.12 

.40 

352 

.92 

.71 

.34 

356 

.90 

• -7.9 

.53 

361 

.92 

.96 

.63 

372 

.91 

, 1.23 

.40 

410 

.91 • 

.81 

.57 

413 

.90 

1.24 

. .32 

424 

.91 

1.61 

.48 

442 

.90 

1.47 

.41 

453 

.91 

.96 

.33 

456 

.91 

.92 

.51 

459 

.91 

1.08 

.31 

471 

.90 

1.43 

.47 

515 

. .90 

1.04 

.25 
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Table 4. 

Spring and Winter Correlations 

of API and Emissivity. 

ft 




Spring 


Winter 



GRID 

Mean 

Emissivity 

Mean 

API 

R 2 

Mean ’ 
Emissivity 

Mean 

API 

R 2 ' 

102. 

.91 

.52 

.22 

.91 

.40 

.04 

106 • 

•>91 

.21 

.10 

.91 

.11 

.05 

133 

.91 

.48 

.06 

.9i 

• .28 

.01 

177 

.91 

.59 

.12 

.91 

.30 

.03 

204 

.91 

.68 

•i 

.24 

.91 

.51 

.11 

213 

.91 

.51 

.15 

.91 

.29 

.04 

240 

.91 

1.08 . 

.37 

.91 

.37 

.04 

256 

.91 

.88 

.08 

.91 

.55 

.26 

259 

.91 

1.15 

. .06 

.90 

.50 

.15 

1 1 288 

.90 

1.08 

.18 

.90 

.73 

..05 

291 

.91 

.89 

.22 

. .91 

.46 

.09 

.'296 

.89 • 

1.29 

.35 

.89 

.44 

.42 

299- 

.90 

1.10 

.11 

.90 

.38 

.19 

343 

.89 

1.34 

.23 

.90 

.57 

.20 

352 

.90 

.95 

.09 

.90 

.63 

.16 

356 

.88 

1.04 

.13 

.88 

.76 

.25 

361 

.90 

1.28 

.13 

.90 

.59 

.09 

372 

.90 

1.27 

.18 

.90 

.54 

.30 

410 ' 

.89 

3.24 

.21 

.89 

-.64 

.17 

413 

.87 

2.09 

.39 

.90 

.73 

.28 

424 

.89 

1.59 

.28 

.89 

.70 

.14 

442 

.87 

2.23 

.43 

.88 

1.15 * 

.31 

453 

.89 

1.20 

.12 

.90 

.89 

.09 

456 

*88 

1.81 

.24 

.89 

.93 

.12 

459 

.88 

2.02 

. ‘.14 

.88 

1.03 

.23 

' 471 

.88 

1.98 

.35 

.89 

.74 

.11 

515 

.88 

3.33 

.38 

.88 

1.71 

.1G 
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The recession coefficient, k, as a function of time was develop 
from a cosine wave to simulate the annual change in evapotrans- 
pi ration. The lowest value of the recession factor was 0.70 in 
August; the highest value was 0.92 in January. 

The exponential 0.891 was used to convert rainfall to effective 
rainfall, which is that portion of rainfall entering the soil for crop 
use. The exponential was applied to the rainfall expressed in mm, 
then the cm equivalent was added to the previous day's index. 

The API model was used for all grid cells without regard to soil 
texture or hydrologic response differences. Soil texture has a 
pronounced effect on microwave response to soil moisture changes 
(Schmugge, 1980; Wang and Schmugge, 1980). With the same moisture 
content by weight, a sandy soil will have a lower emissivity than a 
clay soil due to the greater amount of water in larger pore sizes. 
Schmugge (1980) reported brightness temperature differences of 20 K in 
measurements at 1.55 cm from aircraft sensors over fields in the 
Phoenix AZ area. He reported that conversion of the soil moisture 
parameter to percent of field capacity would normalize each soil 
texture for microwave response to soil moisture. Thus, a geophysical 
data base that includes soil texture would greatly facilitate spatial 
mapping of sol! moisture for crop condition assessment from fallow to 
stand establishment. 

Similarly, the recession factor could be modified to take the 
hydrologic response (e.g., surface runoff and drainage) and the 
evaporation climatology into account for each grid cell. 

To correlate the emissivity with API, four seasons were defined, 
shown in Table 5. The fall (Season 3) correlations are shown in Table 2. 
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Table 5. Seasons tor Emissivity Correlations. 


Season Wheat Development Period 


1 Spring 

2 Summer 

3 Fall 

4 Winter 


major growth 

harvest; fallow 

planting; stand 
establishment 

dormancy 


Feb 1 to Apr 30 
May 1 to Jul 31 

Aug 1 to Oct 31 
Nov 1 to Jan 31 


The number of observations ranged between 64 and 71 for the two years 
(1973 and 1974) of data In the correlations. The coefficients of 
determination (R ) for the test grids range from 0.17 to 0.57, which 
strongly supports a relationship between emissivity and API. The 
significance level is .02 percent for 0.17 and .01 percent for 0.57. 

To determine a sensitivity of the correlations to the size of the 
grid cell, the averaged emissivity for five grids centered upon the 
test grid was correlated with the averaged API. This approximated 
correlations for a 50 km square grid cell, which is the resolution of 
some SMMR and SMMI radiometers. The coefficients of determination 
(R ) presented also in Table 2 do not show any significant differences 
from the R for the grid cell. A very slight improvement was noted 
from most grids. 

The differences in coefficients of determination from one grid 
cell to the next are due to several factors, Including the inade- 
quacies of the API model to describe the moisture content of the emit- 
ting layer. Differences in vegetation that masks the emitting layer 
and hydrological response of each location are primary contributors to 
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the varying correlations. A recession factor for each grid cell 

could be developed to take these factors into account, at least 
empirically. A major conclusion from this aspect of the study is 
that, while spatial mapping of emissivity data will provide qualita- 
tive information on crop moisture, temporal mapping will provide 
quantitative data on crop moisture. Temporal mapping will require the 
refinement of crop moisture models for each general area to take 
vegetation and hydrological response differences into account. 

The summer season coefficients of determination, presented in 

Table 3, show a slightly higher (16 of the 27) level of correlation. 
The number of observations ranged from 38 to 43. 

The winter and spring coefficients are predictably much lower as 
a result of increased vegetation for all sites and the influence of 
snow and frozen ground. These are shown in. Table 4. The number of 
observations in the spring ranged from 69 to 75. The range for winter 
was 82 to 89. 

Snow, ice, and frozen soils have an emissivity similar to that of 

dry soil since the dielectric constant of ice is near that of dry 

soil. In the analysis, the rainfall equivalent of the reported 

snowfall was accumulated then released on the first day that no 
snowcover was reported. The data in Table 4 included all days with a 
reported emissivity. The coefficients of determination improved when 
the analysis was reaccomplished with all days with snow cover 
excluded. For the spring season, the R average increase was 0.10 for 
all 27 test grid cells. Most of the increase was due to a few grid 
cells, which are shown in Table 6. 
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Table 6. Correlations of API and Emissivity for Spring 
with Snow Days Excluded. 


GRID 

R 2 - All days 

R 2 - No snow 

No, .snow days 

213 

.15. 

.31 

21 

• 240 

.37 . 

.58 

22 

356 

.13 

.27 

28 . 

413 

.39 

.57 . 

27 

453 

.12 

.40 

28 

. 456 

.24 

.59 

. 30 

459 

.14 

.33 

31 

. 471 

.35 

.50 

28 

515 

.17 

.38 , • 

21 
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For the winter season, the coefficients of determination for 
several grids showed comparable increases, but some decreases were 
also noted. There was an indication from the air temperature records 
that the snow reports did not represent a uniform snow cover, however. 

A fairly typical emissivity behavior through a period with snow 
and frozen ground is shown in Table 7. Note the rapid and sustained 
increase in brightness temperatures to 0.95 and 0.96 when the tempera- 
tures of the emitting layer were apparently below freezing. The emis- 
sivities remained high without any response to the precipitation that 
was reported on days 365 through 368. The emissivity dropped to 0.86 
after the maximum air temperature climbed well above the freezing 
point. 

Scatterplots for selected grid cells for each season are 
contained in Appendix A. The winter and spring season correlations 
were performed without snow days. 

In a clear, dry atmosphere, the scattering and absorption of 
microwave radiation is negligible. With increasing water vapor, the 
transmissivity decreases, but remains above 90 percent. At 1.55 cm 
wavelengths, a vapor total of two centimeters will attenuate only 
about five percent of the emission (Gloersar and Barath, 1977). Thus, 
for practical purposes, a cloud-free atmosphere is transparent at 1.55 
cm. 

Ice clouds and clouds composed of small water droplets similarly 
are transparent for practical purposes. Larger water droplets in 
clouds and precipitation size droplets (millimeter size) are strong 
reflectors and absorbers of microwave emission, especially in the 
shorter wavelengths. A more-or-less typical thunderstorm with a 
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Table 7. Effect of Low Temperatures on Emissivity 
for a Grid in Harper County, Oklahoma. 



vertical extent of 10 km and a rainfall rate of 5 mm/hr will attenuate 
94 percent of 1.55 cm radiation passing through the cloud. Thus, rain 
clouds In particular will mask the radiation emitted at the surface. 
The liquid water of the clouds and thunderstorm is, however, emitting 
microwave radiation also as a function of temperature and einisslvlty. 
Paris (1969) reported downwelllng radiation at 23 GHz from a heavy 
rain to be 255 k and from a light rain of 200 k, the emissivities of 
rain clouds are of the order of 0.71 to 0.90. These emissivities are 
in the same range as the emissivities of moist soils. 

Thus, thunderstorms in the field of view of ESMR or SMMR will not 
produce a noticeable departure in brightness temperatures unless the 
surrounding area is very dry and hot. The emissivities from a moist 
soil and from a thunderstorm would be very similar. 

In terms of ESMR as an all-weather sensor, the distinction 

between surface emission and atmospheric emission may not be neces- 
sary. The API model based on 24-hour precipitation totals will not 
describe thunderstorms in the field of view, but in general the 
surface and atmospheric emissions are in the same sense for both a dry 
surface - clear atmosphere and, a wet soil - raining atmosphere. 

An examination of several days with extensive, heavy thunder- 
storms at the time of the Nimbus-5 overpass confirms the expected 
effects on the surface microwave emission. On September 26 and 

October 12, 1973, the 17352 radar facsimile chart of the National 

Weather Service showed over five-tenths coverage of moderate or 

greater intensity thunderstorms over south central Kansas and north 

central Oklahoma. The October 12 storms produced especially heavy 

rains, with a 50 cm plus rainfall center near Enid, Oklahoma. The 
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lowest brightness temperature observed tor the flood was 215 k which 
corresponds to an emlssivlty of 0,75 for an emitting temperature of 
287 k. The normalized emlsslvltles for the seven test grid cells in 
the storm area averaged 0,83 on September 26 and 0.82 on October 12, 
Lowest emissivitles were 0.76 on October 12 and 0.81 on September 26, 
Crop Condition Assessment 

The seasonal crop development calendars and the weekly crop 
weather reports for Oklahoma showed normal fall planting and stand 
establishment for the fall of 1974, but extensive delays beginning In 
October 1973. Table 8 contains a summary by week for planting for the 
winter wheat In Oklahoma. 

The delays in seeding, especially in the West Central, North 
Central, and Central crop reporting districts, were also commented 
upon in the weekly crop-weather reports (USOA-USDC Weekly Weather and 
Crop Bulletins). In addition, to the state percentages, the report 
for the week ending October 8, 1973 stated that seeding was 1 to 2 
weeks behind normal. In contrast, the planting activity in the fall 
of 1974 was more-or-less normal, Above normal rains in August 
provided adequate soil moisture for seeding; subsequent rains did not 
cause widespread delays. 

The time series plot of emissivity and API for grid cells 413 and 
343 show the rainfall events that produced the delays. These plots 
are shown in Figures 2 and 3. 

The wet fall that hampered wheat seeding was more pronounced in 
Kansas in 1973. The state average precipitation was 21.21 cm (8.35 
in) compared with 6.71 cm (2.64 in) normal for September. October was 
also wet with 8,15 cm (3.21 in) compared with a normal of 4.75 cm 
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Table 8, 

Seeding Progress 

in Oklahoma for the 

i Fall of 
* 

1973 and 

1974 . 

Crop Reporting 
District 
(and grid Cells) 

Sep 14- 
Sep 15 

Sep 21- 
Sep 22 

Sep 28- 
Sep 29 

Oct 5- ’ 
Oct 6 

Oct 12- 
Oct 13 

Oct 1P- 
Oct 20 

State .Average 
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41 

53" " 

75 


. . 

1974 

11 

31 

38 

65 

88 

95 

Northwest 

1973 

55 

71 

mm mm 

88 

MM 

98 

(213,288,291) 

1974 

41 


93 

99 

MM 

100 

West Central 

1973 

’ . 5 

“ 13 

— , 

52 

MM 

84 

(343) 

1974 

8 

— 

30 

67 . 

mmmm 

. 98 

Southwest 

1973 

6 

11 

MM 

47 

mm mm 

74 

(296,372,424) 

1974 

5 

— 

15 

30 

MM 

86 

North Central 

1973 

5 

15 

MM 

35 

•** rnm 

69 

(413, 515) 

1974 

1 

— 

28 

65 

— \ 

96 

Central 

1973 

7 

24 

MM 

53 

mm mm 

85 

(442, 471) 

1974 

6 

— 

37 

. 63 

MM 

93 
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(1.8/ in). Extracts from the weekly crop weather reports stated: 

September 17 Rains delay planting and seedbed preparation. Planting 
is 10% complete compared with 15% a year ago. 

September 24 Rains delay planting and seedbed preparation. Planting 
is 20% complete compared with 35% a year ago. 

October 1 Considerable reseeding is expected due to heavy rain. 

October 8 Seeding is 10 days behind normal. Seeding is 35% 

complete compared with 85% last year. 

October 15 Seeding is 55% completed compared with 95% last year. 

October 22 Seeding is two weeks behind at 65% complete, compared 

with 100% last year. 

In contrast, 1974 had normal progress statewide through the 
fall. August had 11.10 cm (4.37 in) compared with a normal figure of 
•. 7.65 cm (3.01 in). This guaranteed good planting moisture. September 
was on the dry side, with 4.17 cm (1.64 in) compared with a normal of 
6.71 cm (2.64 in). 

Time series plots for two grid cells in Kansas are shown in 
Figures 4 and 5. In each of the figures, the low emissivities that 
are indicative of the heavy or frequent rains are evident in 1973. 
The contrast with the same period in 19/4 is also evident. A separate 
portion of this research addresses the discrimination capabilities of 
the use of emissivity data to discern precipitation events that can be 
related to progress in field work and stand establishment. 


Microwave Remote Sensing of Crop Stress 

The emitted microwave radiation from the soil is scattered and 
attenuated by vegetation. This is a function of wavelength; the scat- 
tering and attenuation are much more pronounced at the shorter wave- 
lengths such as the 1.55 cm of ESMR and SSMI and the 1.36 and 1.66 cm 
wavelengths of SMMR. As the vegetative cover of the crops increases, 
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the soil component of emission received at a sensor above the canopy 
will decrease. However, the canopy is also emitting radiation as a 
function ot its emissivity and temperature. Mo, et al. (1981) found 
the effective canopy thickness to be directly proportional to the 
amount of water present in the plant materials. 

For a crop with adequate moisture, the canopy temperature will be 
at or below the air temperature in mid-day due to the cooling effects 
of evapotranspiration. Further, the high moisture content of the crop 
will decrease the emissivity. Thus, the brightness temperature of a 
well-watered crop will be lower than that of a crop experiencing 
moisture stress. In addition, plants under moisture stress often 
•« exhibit leaf rolling or wilting, which decreases the attenuation and 
scattering. More of the soil surface is also exposed, so the net 
effect is for a greater component. of emission from the dry soil to 
reach the sensor. The net effect is for the brightness temperatures 
to respond in the same sense for a well-watered crop and an moist soil 
surface. Low brightness temperatures and low emissivities are indica- 
tive of adequate crop moisture, while high brightness temperatures and 
high emissivities are indicative of either a dry soil surface that, if 
persistent, could indicate insufficient crop moisture for plant and 
stand establishment or a crop canopy that is experiencing moisture 
stress. 

The emissivities for grid cells with large winter wheat acreages 
conformed fairly well with expectations. For the months of April and 
May, the winter wheat canopies will be at their maximum extent. For a 
vigorous crop without moisture stress, the emissivity should be fairly 
constant with a very strong correlation with soil moisture, as 
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inferred by API. Table 9 shows the comparison of two grids in the 
winter wheat areas of Oklahoma for 1974 and 19/5. Two periods were 
selected with varying rainfall and presumably soil moisture during 
periods when the winter wheat canopy should be fully developed. 

For grid 413 with apparently adequate soil moisture, the 

emissivities were fairly constant in 1974 and 1975. The emissivities 
for corresponding periods for grid 213, with much less rainfall, were 
higher. Indications of moisture stress at these times are shown by 
the yields. 

The correlation coefficients between API and emissivity for these 
grid cells for the periods show a range from -0.12 to -0.64, which is 
not surprising in view of the number of factors involved. The 

response of emissivity to rainfall, as indicated by the greater 

negative correlation coefficients, will be from several sources. Soil 
with little or no vegetative cover, emission and scattering in the 
atmosphere, and limitations of the simple models used will all 

contri bute. 

The major conclusions of this aspect of the investigation are: 

1. A we II -watered crop will have a decreased emitting 
temperature due to ET and a reduced emissivity due to an increased 
water content of the canopy. Both contribute in the same direction 
toward a decreased brightness temperature when compared to a crop 
under stress. The ESMR data support this concept. 

2. The model developed by Blanchard, _et _al_. (1981) has the 
potential for use as an early warning for moisture stress. The 
summation of daily departures of emissivity from the well-watered 
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canopy emissivity (of 0.90) will be directly proportional to 
accumulated moisture stress of the crop. This should serve as an 
early warning screening device for yield reduction. 

3. The simple models used have significant utility for an 
all-weather crop condition screening device for large relatively 
monoculture agricultural areas. 

Table. 9 ESMR Microwave Emissivities from 
Vegetated Terrain. 


Grid 

(Jounty(OK) 

Planted Acreage 

Yield 



(1000s ha) 

(T/ha) 



1973-74 1974-75 

'1973-74 

213 

Beaver 

122 131 

0.6 

413 

Alfalfa 

125 133 

1.9 


1974-75 


0.9 

1.9 


April 11 - May 6, 1974 May 4 - May 30, 1975 

Grid 213 


Number of Emissivity Obs. 15 15 

Emissivity mean .934 .916 

Standard deviation .020 .026 

Total rainfall (cm) 0.10 9.45 

API - emissivity correlation -0.12 -0.64 


Number of Emissivity Obs. 
Emissivity mean 
Standard deviation 
Total rainfall (cm) 

API - emissivity correlation 


Grid 413 
15 

.896 

.011 

12.81 


-0.29 


16 

.897 

.019 

22.02 (Apr. 

27 -May 30) 
-0.51 
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RELATIVE FREQUENCY DISTRIBUTIONS 

To provide a statistical approach to the identification of dry 
periods, a number of cumulative relative frequency distributions were 
computed. Four are shown in Figures 6 through 9. These frequency 
distributions relate emissivity to the precipitation history in a 
probabilistic manner. 

The distributions were constructed for grid cell data using all 
27 case study grid cells for season 3 (Aug, Sep, Oct). The following 
variables were defined. 

Variable 
RAINO 
RAIN5 
RAIN10 
RAIN15 
RAIN20 

Additionally, five emissivity variables were defined. 

Variable 
ECATO 
ECAT5 

ECAT10 

ECAT15 

ECAT20 

These emissivity categories were defined, to be determined from the 
minimum emissivity for the period defined: 


Definition 

Emissivity category for day of observation 

Lowest emissivity category for day of 

observation and previous 5 days 

Lowest emissivity category for .day of 

observation and previous 10 days 

Lowest emissivity category for day of 

observation and previous 15 days 

Lowest emissivity category for day of 

observation and previous 20 days. 


Definition . 

Precipitation on Day of Observation 
RainO + Amount for Previous 5 days 
RainO + Amount for Previous 10 days 
RainO + Amount for Previous 15 days 
RainO + Amount for Previous 20 days 



... 


28 


ORIGINAL PAGE m 
OF POOR QUALITY 


Emissivity Category 
1 
2 

3 

4 

5 

6 

7 

8 


Definition 
min(EMIS) < 0,78 
0.78 < min(EMIS) _< 0.81 
0.81 < min(EMIS) 1 0.84 
0.84 < min(EMIS) .< 0.87 
0.87 < min(EMIS) < 0.90 
0.90 < min(EMIS) < 0.93 
0.93 < min(EMIS) j< 0.96 
0.96 < min(EMIS) 


Each of Figures 6 through 9 have several cumulative distribution 
curves, one for each of several values of ECAT. The interpretation of 
these curves is straightforward. For instance, Figure 6, the 
distribution of RAIN5 by ECAT5, shows that 90% of the days with 
ECAT5=5 (0.87 < min(EMIS) £ 0.90) had six=day rainfall of 1.3 cm or 
less. If these relative frequency distributions can be assumed to 
approximate true probabilities, a 6-day minimum emissivity higher than 
0.90 (ECAT0 = 6) would indicate a probability of 99% that the rain for 
the day of observation and previous five days was no more than 1 cm. 

The spacing of curves on these graphs indicates that emissivity 
data can be used to discriminate between relatively moist and dry crop 
soil moisture conditions; if this were not the case, all of the curves 
on any of the figures would approximately coincide with all of the 
others. In fact, the curves are nearly parallel and have 

approximately regular spacing over most of the emissivity range. It 
is only for the extremes of emissivity, where the sample sizes are 
small, that these relationships are not valid. Apparently, the true 
probability curves approximated by these sample distributions are 
well-behaved smooth curves arranged in order by the value of ECAT. 
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Figure 6. Cumulative Frequency Distribution of RAIN5 by ECAT5 for Season 3 (Aug, Sept, (Oct) 
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Figure 7. Cumulative Frequency Distribution of RAIN10 by ECAT10 for Season 3 (Aug, Sept, (Oct). 
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This behavior extends throughout the series of figures, giving 
the capability to detect extended periods of dry conditions. For 
instance, Figure 9, the distribution of RAIN20 by ECAT20, indicates 
that a day with ECAT20 = 5 (min 21-day emissivity between 0.87 and 
0.90) has a 90% probability of being associated with a 21-day rainfall 
of no more than 6 cm. 

The statistical behavior expressed in these curves can be under- 
stood from the physics of soil moisture. As shown in the time series 
plots, Figures 2 through 5, each rainfall event diminishes the emis- 
sivity of the emitting layer by an amount proportional to the effec- 
tive precipitation. The soil then undergoes relatively slow stage 1 
drying for summer days with a correspondingly slow emissivity 
increase, followed by faster stage 2 drying with rapid return to the 
original emissivity. The length of stage 1 drying makes these changes 
detectable, even with ESMR's irregular observation frequency. In 
general, then, the emissivity alternates between a dry surface condi- 
tion at about 0.92-0.94, and inoister conditions with lower enissivi- 
ties indicative of rainfall amount. Superimposed on this behavior is 
a random noise component resulting from irregularities in observation 
frequency, soil type, crop development stage, system noise, 
unaccounted meteorological effects,' etc. Even with this high-ampli- 
tude noise, if the lowest emissivity in a given period is sufficiently 
high, it can be concluded that rainfall during the period was very 
low. With decreasing minimum emissivity, the estimate of total rain- 
fall must increase. 

These relationships have been summarized in Figure 10, which 
shows the 75% confidence level for precipitation in any of several 
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0 1 2 3 4 5 6 7 

ECAT 

A RAIN 5, GIVEN ECAT 5 

B RAIN 10, GIVEN ECAT 10 

C RAIN 1 5, GIVEN ECAT 15 
D RAIN 20, GIVEN ECAT 20 


Figure 10. 75% Confidence Level for Total Rain in Period, Given 

Emissivity Category (ECAT) of Lowest Emissi vity . 
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periods, given the minimum emissivity of the period. For instance, a 
minimum emissivity of 0.88 for any 21-day period (ECAT20 = 5) 
indicates, with 75% confidence, that rainfall (RAIN2Q) in the same 
period was no more than 4.7 cm. Figure 10 can, thus, be used to 

identify periods of low rainfall. 
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CLASSIFICATION BY DISCRIMINANT ANALYSIS 


Passive microwave data were used in discriminant models to 
classify the moisture conditions in wheat croplands. Three elements 
are required: an objective set of moisture categories defined by 

ground, truth data, one or more remotely sensed variables 
characteristic of each moisture category, and a method of analysis by 
which to relate the two. The method of discriminant analysis was used 
successfully. 

Moisture Categories 

Two methods of categorization were employed. In the first. 


categories were defined by API: 


API CAT Definition 

1 API ,<1.0 cm 

2 1.0 cm < API j, 2.5 cm 

3 2.5 cm < API 


This system was used only for the months of March, April and May. 

In the second system, moisture categories were defined by total 
precipitation amounts received over the three-week period ending on 
the day Ov the observation. To examine the effect of the number of 
categories used, three levels of detail were analyzed. 


RAINCAT1 

1 

2 

3 

4 


Definition 
Rain = 0 cm 

0 cm < Rain £ 1 cm 

1 cm < Rain _< 5 cm 
5 cm < Rain 
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RAINCAT2 

1 

2 

0 cm 

3 

1 cm 

4 

5 cm 

5 

15 cm 

RAINC/T3 


A 

2 

0 cm 

3 

1 cm 

4 

3 cm 

5 

5 cm 

6 

TO cm 

7 

15 cm 

8 

20 cm 


Definition 
Rain = 0 cm 

< Rain < 1 cm 

< Rain < 5 cm 

< Rain < 15 cin 

< Rain 

Definition 

Rain = 0 cm 

< Rain £ 1 cm 

< Rain _< 3 cm 

< Rain £ 5 cm 

< Ra'i n <_ 10 cm 

< Rain j< 15 cm r > 

< Rain <. 20 cm 

< Rain 
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These RAINCAT systems were used for each of the four seasons defined 
as: 


SEASON 

1 

2 

3 

4 


MONTHS 

Feb-Mar-Apr 

May-Jun-Jul 

Aug-Sep-Oct 

Nov-Dec-Jan 


Classification Variables 


The following variables were used: 

VARIABLE DEFINITION 

EMI SO Emissivity on day of observation 

EMIS5 Average emissivity for 1st through 5th days 

before observation. 

EMIS1Q Average emissivity for 6th through 10th days 

before observation. 

EMIS15 Average emissivity for 11th through 15th days 

before observation. 

EMIS20 Average emissivity for 16th through 20th days 

before observation. 

EMISO was used as the single classifying variable in the APICAT 

# 

study. All five variables were used in the RAINCAT classifications. 
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Two additional analyses were performed, to evaluate the utility 
of restricted variable sets. These analyses were motivated by the 
method's handling of missing data for the RAINCAT3 analysis. For 
example, a 5-day break in the ESMR emissivity data might result in, 
say, EMIS 5 being undefined while the other four variables were 
available for anlaysis. In this case, the discriminant analysis would 
simply ignore the entire observation vector. On the other hand, a 
separate discriminant analysis could be performed, using all of the 
variables except EMIS5. To examine the effect of restricted variable 
sets, the following analyses were performed, using the same categories 
as in RAINCAT3: 

RAINCAT4 Variables: EMISO, EMIS5, EMIS10 

RAINCAT5 Variables: EMISO, EMIS10, EMIS15, EMIS20 

Analysis Method 

The theory of discriminant analysis is presented in such texts as 
Kendall and Stuart (1976) and Rao (1965). The method has been 
implemented in the Statistical Analysis System (SAS) as Procedure 
DISCRIM (SAS Institute Inc., 1979). SAS is a software package 
available on the AMDAHL 470 computer system. Inputs to the DISCRIM 
procedure are observations of classification variables and related 
true categories determined from ground truth data. The DISCRIM 
procedure in SAS develops a system of probabilistic discriminant 
functions. It applies these functions against the input data set and 
prepares a summary of their performance in classifying observations 
into categories. Optionally, SAS can store these functions for later 
use in classifying independent data sets. 
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DISCRIM works by assuming that the observation vectors from each 
category are samples from distinct multivariate normal populations. 
For cateogry t of n categories, it calculates a mean n-dimensional 
observation vector It then computes a set of generalized squared 

distance functions, D£(x), characterizing the separation in n-space of 

each vector x from’Xj.. The form of the distance function depends upon 
the homogeneity of the m within-groups covariance matrices, and upon 
the assumed prior probabilities. If the several within-category 

populations can be assumed to share a common covariance matrix, it is 
estimated by a pooled covariance matrix S, computed from all of the 
observations. If this assumption is not made, separate within-groups 
covariance matrices, St, are computed. DISCRIM also has an optional 
test of the covariance homogeneity hypothesis. The distance function 
has the linear form: 

d|(>c) = g x (x,t) + g 2 (x,t) 

where 

g : (x,t) = (x - x t )'Sj 1 (x-x t ) + ln|s t | 

if within group covariance matrices are used. 

g^x.t) = (x - x’ r )'S“ 1 (x - ;x t ) otherwise. 

g 2 (x,t) = -2 In (prior probability for group t), 

if prior probabilities are not assumed equal. 
g 2 (x,t) = 0 otherwise. 
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The classification procedure, then, is to assign each observation 
to group u if D£(x) is a min ’mum for t=u. DISCRIM takes this process 
one step further, by using the assumption of multivariate normality to 
compute posterior probabilities for membership in each of the in 
categories. These probabilities are given by 

P(j|x) = exp[-|o?(x)]/ l exp(-7rD^(x) J 

k “■ 1 

From this probability information, a user can decide how much 
confidence to place in the system's classification of an observation. 
Table 9 is an example of DISCRIM's output listing of classifications 
and posterior probabilities for a calibration data set. Errors in 
classification are flagged by an asterisk. 

Prior Probabilities 

The discriminant method requires estimation of the prior 
probability of each category; that is, the assumed probability that an 
event will fall in a given category, based on all that is known prior 
to acquisition of the values of the classification variables. 
I) I SCRIM can use any set of assumed prior probabilities, but two such 
sets were used. One results from the assumption that, until current 
remotely sensed input data are available, the probability of each 
category is equal to that of any other category. This option will be 
referred to as "equal priors". The other, perhaps more reasonable 
assumption would be that prior probabilities are equal to the observed 
relative frequencies of the categories. This option is referred to as 
"proportional priors". 
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Table ID Example of Output from DISCRIM Procedure 


PROPORTIONAL PRIOR PRORAkUIT/lS 
AVlO 

CLASSIFICATION RESULTS POR CAL I RRAT ION DATA 


POSTERIOR PRORAR I L I TV OP MSMSERSNIP IN MO| STCAT $ 


to 

PROM CLASSIPISO 

1 

3 

3 



1401 STCAT INTO MOISTCAT 




f ! 

31(741 

1 

1 

0,4471 

0,1411 

0,0014 

S! 

H 

311743 

1 

1 

0,1333 

0,3334 

0,0130 

a 

111744 

1 

I 

0,4041 

O, 1744 

0,0131 

* 

111741 

3 

1 • 

0, 111 1 

0,1173 

0,0431 

j| 

USTlO 

3 

3 

0,0131 

0.3131 

0,1033 


RSRSRS 

3 

1 • 

0,4701 

0,4304 

0,1031 

3M2I7 

a 

1 • 

0.4144 

0,1317 

0,00*1 

| 

311311 


1 

0,4000 

0, 1411 

0,0134 

1 

311 300 

1 

1 

0,4144 

0,4044 

0,0110 

311301 

1 

1 

0,7771 

0,3044 

0.0173 

i 

311303 

1 

1 

0* 4373 

0.3137 

0,0410 

! 

311304 

3 

3 

0,0314 

0,4341 

0,1341 

f! 

311301 

3 

3 

0,0171 

0,1314 

0,3773 

ii 

311301 

3 

3 

0, 104 | 

0,1344 

0,3171 

)! 

311311 

1 

a • 

0.3114 

0.4401 

0. 114 1 

ri 

311313 

1 

i 

0,7413 

0.3301 

0,0333 

ii 

312)3 14 

1 

i 

0,4443 

0,8714 

0,0311 

■j 

311333 

1 

i 

0,7334 

0,3444 

0.0371 

ii 

311333 

1 

i 

0,7441 

0,3344 

0,0337 


311331 

1 

i 

0,7137 

0. 1114 

0,0141 


311331 

1 

i 

0,4143 

0.1347 

0,0013 

* 

311340 

1 

i 

0,1114 

0.3133 

0,0813 


311343 

1 

i 

0.7031 

0,31«1 

0,0331 


311344 

1 

i 

0,7317 

0.3431 

— ,£>253 


311341 

1 

i 

0,1737 

0,3441 

0,031* 


311341 

1 

i 

0,4113 

0, 1301 

0,0031 


311347 

1 

i 

0,4033 

0,1*47 

0,0131 


311341 

1 


0,4114 

0. 1311 

0,0047 


31131! 

1 

i 

0,4717 

0.3483 

0,0400 


31*313 

1 

i 

0,4317 

0. 1111 

0,0041 


311313 

1 

i 

0,4141 

0, 1370 

0,0041 


31SM4 

» 

i 

0,4341 

0* 114 | 

0,0013 


311 L ^1 

1 

i 

0.4143 

0. 1314 

0,0033 


311314 

1 

i 

0,4131 

O, 1741 

0,0113 


311311 

1 

i 

0,4411 

0. 1444 

0.00* 1 


311311 


i 

0,7117 

0. 1131 

0,0148 


311313 

1 

i 

0,4390 

0, |130 

0,0010 

£ 

311311 

1 

i 

0, 404 1 

0.3341 

0,0174 

* ;■ 

311311 

1 

i 

0.4114 

0. 1313 

0,0013 


311314 

1 

i 

0,4414 

0.1113 

0,0073 


31137 1 

1 

i 

0.4443 

0, 1341 

0,003 1 


311373 

t 

i 

0,4343 

0, 1171 

0,0083 


311373 

1 

i 

0,7474 

0,3133 

O.Otl 1 


311371 

1 

i 

0, 4444 

O, 1403 

0,0014 


311377 

3 

! • 

0,4140 

0,1731 

0.0111 


311374 

3 

1 • 

0»£041 

0, 1137 

0*0137 


311371 

3 

1 • 

0.4141 

0, 13*1 

0,0034 


311343 

3 

1 • 

0,7144 

0,3134 

0.0380 


311344 

3 

1 « 

0,4131 

0,1741 

0.0114 


311341 

1 

1 

0,4314 

0,111* 

0.0088 
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Procedure 

Analyses were performed for grid cell data representing 25 km 
resolution elements and for spatially averaged data representing 
simulated resolution of about 50 km. These data types are denoted, 
respectively, by "AV = 0" and "AV = 1". 

The APICAT analysis was performed for each of the four combi- 
nations of averaged/unaveraged data and equal /proportional prior prob- 
abilities; and the period of study was the months of March through May 
of 1974 and 1975. The RAINCAT analyses were performed for all com- 
binations of season number and equal /proportional priors, but only the 
non-averaged data set was used. In each case, days with snowcover 
were allowed to remain in the data. 

In every analysis, the covariance homogeneity test rejected (at 
a = 0.10) the hypothesis of homogeneity. Thus, individual 

within-groups covariance matrices were computed. 

Results 

Three measures of effectiveness were used to evaluate the various 
analyses. Prefigurance (PF) and Post Agreement (PA) (AWS, 1978) mea- 
sure the conditional probability of a correct classification given, 
respectively, the observed category or the classified category. For 
instance, suppose category 1 occurs 10 times, 20 observations are 
classified as category 1, and only 5 of those classifications are 
correct. The prefigurance for category 1 is 5/10 = 50%. The post 
agreement for category 1 is 5/20 = 25%. Each category, then, has both 
a PF and a PA score. Appendices C and D give prefigurance and post 

agreement results for each analysis in a relative frequency matrix. 
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The individual category PF or PA scores are the numbers on the 
diagonals of the matrices. An overall PA or PF score can be defined 
for each matrix by taking the simple average of all defined scores for 
that matrix. For example, if a 5 x 5 matrix has one undefined 
diagonal entry, and four defined PF scores, then 

l 4 

W » t I PF i 
4 i=l 1 

This average score must be used with caution, since it gives undue 
weight to Individual scores from categories with low frequencies of 
occurrence. Nevertheless, mean prefigurar$Te and postagreement scores 
can be used to compare the performance of different models. 

The third measure of classification effectiveness is simply the 
relative frequency of correct classifications, denoted "%C0R". This 
measure and both "PF and ¥K are given in Table 10 for each analysis. 

Discussion 

These summary figures for the APICAT analyses show no appreciable 
difference betwen results for averaged and nonaveraged data. This 
result will be important to the SMMR phase of this project, as it 
demonstrates that reduced resolutions have little effect on 
classification accuracy. 

A second observation is that, for all analyses, the assumption of 
proportional prior probabilities results in better total percent 
correct (%C0R) as well as better mean post agreement (PA). By con- 
trast the assumption of equal priors gave equal or, usually, better 
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Table 11. Performance Summary 


API CAT (Mar-Apr-May) 

%C0R 

PF 

PA 

N 

AV = 0, PE* 

62 

58 

54 

2165 

AV » 0, PP** 

65 

53 

60 

2165 

AV = 1, PE 

62 

58 

53 

2228 

AV = 1, PP 

64 

51 

58 

2228 

RAINCAT1 (AV = 0) 

%C0R 

PF 

PA 

N 

Season 1, PE 

36 

46 

41 

1308 

PP 

59 

39 

53 

1308 

Season 2 , PE 

64 

65 

59 

769 

PP 

74 

65 

70 

769 

Season 3, PE 

65 

76 

64 

1330 

PP 

70 

76 

69 

1330 

Season 4, PE 

47 

58 

48 

1818 

PP 

55 

47 

56 

1818 


RAINCAT2 (AV = 0) 

%C0R 

PF 

PA 

N 

Season 

1, 

PE 

32 

52 

34 

1308 



PP 

59 

36 

54 

1308 

Season 

2, 

PE 

56 

68 

48 

769 



PP 

72 

54 

77 

769 

Season 

3, 

PE 

53 

69 

55 

1330 



PP 

63 

67 

67 

1330 

Season 

4, 

PE 

44 

59 

39 

1818 



PP 

54 

39 

52 

1818 


Continued 


*PE - Priors equal 
**PP - Priors proportional 
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fcCOR PF PA 

N 

Season 1, PE 

29 

36 

29 

1308 

PP 

49 

22 

42 

1308 

Season 2 , PE 

39 

44 

40 

769 

PP 

55 

35 

55 

769 

Season 3, PE 

36 

50 

39 

1330 

PP 

47 

48 

50 

1330 

Season 4, PE 

34 

58 

37 

1818 

PP 

48 

40 

55 

1818 

RAINCAT4 (AV = 0) 

%C0R 

PF 


N 

Season 1, PE 

27 

41 

22 

1645 

PP 

42 

27 

39 

1645 

Season 2, PE 

32 

38 

34 

928 

PP 

49 

25 

45 

928 

Season 3, PE 

25 

41 

25 

1629 

PP 

40 

28 

35 

1629 

Season 4, PE 

29 

48 

29 

2052 

PP 

48 

25 

43 

2052 


RAINCAT5 (AV = 0) 

%C0R 

PF 

PA 

N 

Season 

1. 

PE 

24 

55 

30 

1362 



PP 

47 

32 

63 

1362 

Season 

2, 

PE 

32 

38 

33 

769 



PP 

52 

29 

58 

769 

Season 

3, 

PE 

33 

48 

33 

1357 



PP 

45 

38 

50 

1357 

Season 

4, 

PE 

32 

55 

32 

1899 



PP 

46 

33 

47 

1899 
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SUMMARY 

This research was conducted to examine the potential of the use 
of short wavelength information from passive microwave radiometers on 
earth satellites. The Electrically Scanning Microwave Radiometer 
(ESMR) 1.55 cm passive microwave radiation, expressed as a brightness 
temperature, was converted to emissivity for 25 cm grid cells for 
September 5, 1973 to May 30, 1975 for the Southern Great Plains. The 
frequency of coverage was on the order of once every two or three days 
for the majority of the period for the eastern two-thirds of Kansas 
and Oklahoma and northwest Texas. Daily estimates of air temperature, 
precipitation, and snow cover ware also available for the grid cells. 
Correlations of these emissi vitles showed the following results: 

1. ESMR emissivities were highly correlated with an antecedent 
precipitation index (API) used to infer the moisture content of the 
upper layer of the soil. Correlations were highest in the grid cells 
with high percentages of winter wheat in the fall at planting time. 

2. Temporal series of ESMR emissivity related well with crop 
calendar documentions of the progress of planting, the state of soil 
moisture at planting, and the occurrence of excessive moisture that 
necessitated replanting. 

3. Case study analyses of emissivity identified periods of 
frozen soils, as inferred from air temperature records. The possibil- 
ity of winter kill detection or early warning is suggested. 

4. Emissivities for fully developed winter wheat crop canopies 
in April and May correlated well with adequacy of crop moisture, as 
indicated by rainfall reports and yield tabulations. Crop canopies in 
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a fully-watered state were lower (0,897) than canopies under stress 
(0.930). Small standard deviations Indicate the significance of the 
small variations. 

5. The emlssivitles from thunderstorms In the grid area at the 
time of overpass were of the same range as emlsslvlties for wet soil. 
Thus, the presence of thunderstorms will not significantly degrade the 
ability to monitor crop moisture with passive microwave data. 

6. Observed emissivities from ESMR were very similar to the 
emi ssivities obtained form aircraft and truck measurements. This 
indicates the validity of approximating the temperature of the emit- 
ting layer with surface reports of air temperature. 

7. Passive microwave sensors from space nave all-weather, day or 
night utility. 

8. The use of 1.55 cm passive microwave emissivities for early 
warning of crop condition assessment is strongly supported by the 
investigators. 
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APPENDIX A 

SCATTER PLOTS OF EMISSIVITY AND API FOR 
SELECTED GRID CELLS FOR EACH SEASON 


STATISTICAL ANALYSIS SYSTEM 

SERS0N*2 GHI D=2U0 




SYSTEM 


T A T I 


S T 1 C A L A N A L Y S 

SERSQN»3 GniD“240 
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s 



S y S T E M 


STATISTICAL ANALYSIS 

S£flSQM»U GPID«2*iO 
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STATISTICAL ANALYSIS SYSTEM 

SERSOH-1 G R 1D“413 
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STATISTICAL ANALYSIS SYSTEM 

SEflS0N=2 GRID»413 
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STATISTICAL ANALYSIS SYSTEM 

SEfiS0K«3 GRin=m3 
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S T A T I S 


TIC AL ANALYS 

SEASON GRI D°4 1 3 


S SYSTEM 
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STATISTICAL ANALYSIS SYSTEM 

SERSGN“ 1 GBID“45S 



e?i 
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STATISTICAL analysis system 

SERSOf.'nS GRIO^SB 
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TATISTICAL ANALYSIS 

SF.RS3N-3 CRID»456 
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SYSTEM 





STATISTICAL ANALYSIS SYSTEM 

SEBSON-M GRID=y55 
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STATISTICAL ANALYSTS SYSTEM 

SEfiS0N«= 1 GR1 0=51 5 


ed: PS 
1.0 7S-J 


i.gso- 


1.035- 


I. 000- 
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0.9754 



m 
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STATISTICAL ANALYSIS SYSTEM 

SEASON-2 r -,\ '515 • 
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STATISTICAL A N A L Y S I S 

SERS0U=3 GHID=515 



0. 77S«] 


S Y STE M 


0* 7S04J 


1 2 $ 


<3 


wyiu . - iv^wyav., acyywygr 

S 6 7 o 9 

fiP! 


.v.\ t vygg vrr in mwi yywi^-s - s ’ ggy 

10 II 12 13! 


66 


STATISTIC A L ANALYSIS S Y S T E M 


SEASONS GRID-515 



%»*»** 
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APPENDIX B 

DISTRICT PRECIPITATION AVERAGES (CM) FOR 
KANSAS AND OKLAHOMA FOR JUNE TO NOVEMBER 1973 AND 1974 
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Month District 



November 


ndftlNAL 


APPENDIX C 

PREFIGURANCE TABLES* 


*Relative frequency of actual category (from ground truth), given the 
category into which observations were classified. If the actual cate 
gory is a and the classified category is c, these tables give 
approximations of the conditional probabilities, P(c a). Prefigur- 
ance, for any category, is the entry on the diagonal. 
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KA INCAT 1 





2 

3 

4 

N 

23 

6 

3 

66 

36 

7 

6 

486 

17 

27 

31 

593 ii 

ti 

11 

21 

53 

165 | 



1308 j 



2 

3 

. 4 

N 

85 

16 

2 

66 

CO 

CO 

16 

1 

486 

37 

59 

4 

593 J 

,j 

26 

1 62 

12 

163 j 

! f 

, i 


1308 j 
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PREFIGURANCE TABLES 


RAINCAT 1 


PRIORS EQUAL 



1 

2 

3 

4 

1 

- 

- 

- 

- 

2 

0 

76 

22 

2 

3 

* 

0 

41 

40 

• 19 

4 

0 

3 

17 

80 


* 


N 

0 

90 

290 

389 

769 


PRIORS PROPORTIONAL 



1 

2 

3 

4 

1 

0 

0 

0 

o 

2 

0 

38 

59 

3 

3 

0 

7 

73 

20 

4 

0 

*0 

16 

84 


N 

0 

90 

290 

389 

769 


SEASON = 2 


73 
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PREFIGURANCE TABLES 


R A INCAT 1 
PRIORS EQUAL 



1 

2 

3 

4 

1 

100 

0 

0 

0- 

2 

2 

82 

15 

2 

3 ' 

o' 

31 

52 

17 

4 

0 

4 

27 

69 


N 

8 

168 

448 

706 

1330 


PRIORS PROPORTIONAL 



1 

2 

3 

4 

1 

100 

0 

0 , 

0 

2 

1 

" 71 

24 

3 

3 

0 

19 

58 

23 

4 

0 

1 

23 

76 


N 

8 

168 

448 

706 

1330 


SEASON » 3 


74 



RA INCAT 1 
PRIORS EQUAL 


PREF1GURANCE TABLES 



1 

2 

3 

4 

1 

76 

9 

13 

3 

2 

38 

27 

28 

6 

3 

15 

11 

55 

19 

4 

2 

2 

21 

75 


N 

150 

761 

728 

179 

1818 


PRIORS PROPORTIONAL 



1 

2 

3 

4 

1 

21 

40 

39 

0 

2 

7 

41 

51 

1 

3 

1 

16 

79 

3 

4 

0 

3 

51 

46 


N 

150 

761 

728 

179 
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PREFIGURANCE TABLES 

RA INCAT 2 
PRIORS EQUAL 



1 

2 

3 

4 

5 

1 

68 

18 

6 

3 

8 

2 

51 

32 

6 

7 

4 

3 

25 

16 

22 

29 

9 

4 

15 

10 

17 

52 

5 

5 

0 

13 

0 

0 

88 


N 

66 

486 

593 

155 

8 

1308 


PRIORS PROPORTIONAL 



85 

14 

2 

0 

0 

83 

16 

1 

0 

0 

37 

59 

3 

«0 
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naiGINM" PA QE JS 


APPENDIX D 

POSTAGREEMENT TABLES* 


*Relati ve frequency of classification in each category, given the 
category into which observations were classified. If the actual ca 
gory is a and the classified category is c, these tables 9™®. 
approximations of the conditional probabilities, P(a|c). Postagree- 
ment, for any category, is the entry on the diagonal. 
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